The Effect of Laser Focusing Conditions on Propagation and Monoenergetic Electron 

Production in Laser Wakefield Accelerators 
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The effect of laser focusing conditions on the evolution of relativistic plasma waves in laser wake- 
held accelerators is studied both experimentally and with particle-in-cell simulations. For short focal 
length (wo < A p ) interactions, beam break-up prevents stable propagation of the pulse. High held 
gradients lead to non-localized phase injection of electrons, and thus broad energy spread beams. 
However for long focal length geometries (u>o > A p ), a single optical filament can capture the ma- 
jority of the laser energy, and self-guide over distances comparable to the dephasing length, even 
for these short-pulses (cr m Xp). This allows the wakefield to evolve to the correct shape for the 
production of the monoenergetic electron bunches, as measured in the experiment. 
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The concept of using lasers to accelerate particles to 
relativistic energies in plasma [l| has been brought closer 
to practical realization by recent experimental results. 
These have demonstrated the production of relativis- 
tic electron beams with high energies @, 0] , low energy 
spread [3, 0, 0] and low emittance Q. These accelera- 
tors promise to revolutionise the many uses of such par- 
ticle beams, due to their high accelerating fields and thus 
compact size. Of the numerous laser based schemes for 
plasma accelerators [l|, 0], only the laser wakefield ac- 
celerator (LWFA) scheme has been capable of producing 
quasi-monoenergetic electron beams [I0,0]- The LWFA 
consists of a high intensity laser propagating through a 
plasma, such that the pulse length is shorter than the 
relativistic plasma wavelength (ct < X p = 2ttc/u) p ). 

The longitudinal electric field of the wakefield — a rel- 
ativistic electron plasma wave, with phase velocity (3 p hC, 
generated by the ponderomotive force of the pulse — can 
be used to accelerate trapped electrons. The maximum 
electric field the wave can support, in the cold one- 
dimensional limit, is E max = [2(7ph — 1)] 1 ' 2 -Eq) where 
E a = m e cu; p /e and j ph = (1 - P ph )~ 1/2 M- when the 
wave reaches this wavebreaking threshold, the maximum 
longitudinal velocity of plasma wave electrons exceeds 
PphC. Electrons are self-trapped and accelerated forward, 
instead of oscillating around their initial position. 

In multiple dimensions the self-trapping mechanism 
and the wakefield dynamics can be very different from the 
one-dimensional case [ToL 11, 12\. In three dimensions, 



due to the large radial component of the electron trajec- 
tories, wake oscillations are heavily damped [l3j to the 
extent that in the extreme case there is only one plasma 
wave period 14[. Due to electron cavitation and ion in- 



ertia, the electric fields in this wave period are approxi- 
mately linear and focusing towards the centre of the ion 



cavity. Electrons slipping back with respect to the pulse 
are injected into the rear of the cavity where the electric 
field is strongly forward accelerating. If the charge in the 
injected bunch is high enough, its own electric field coun- 
teracts that of the cavitated region, preventing further 
injection and creating a short, monoenergetic bunch. 

The electrons in the bunch can reach energies deter- 
mined by the electric field in the plasma wave and the 
length over which the acceleration takes place. The max- 
imum energy that an electron can gain is set by the de- 
phasing length Ld ~ 7 p h^pi which is the distance before 
the relativistic electrons out-run the accelerating phase of 
the wake [H . For this to happen the laser must propagate 
at high intensity to the order of this dephasing length, 
which is generally longer than its Rayleigh length, zr. 

Due to the increased energy requirement and technical 
difficulty in producing intense, ultra-short pulses, initial 
LWFA studies used short focal length optics to obtain in- 
tensities sufficient to reach the wavebreaking threshold. 
(A criteria for injection is that the laser normalized vector 
potential ao = eA/m e c > 4.3 [13]). However, experimen- 
tal results with tight focusing have been less successful in 
terms of charge, electron energy and spectral shape than 
with longer focal length focusing optics and ao ~ 1 [3, 0] ■ 

In this letter we examine the importance of the fo- 
cusing geometry on the evolution of plasma waves and 
laser pulse propagation in a LWF accelerator. This is 
done experimentally by directly comparing laser pulse 
and electron beam characteristics using two different fo- 
cusing geometries, chosen to produce either a spot of size 
Wq < A p or wo > X p for ct < X p . We demonstrate that 
for wq > A p , self- focusing permits ao to be sufficiently 
amplified to allow self-injection. Moreover the quasi- 
static nature of this compression allows for an adiabatic 
(i.e. controlled) approach to self-injection, which results 
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in production of monoenergetic bunches. We also show 
that the pulse is maintained at high intensity, through 
self-guiding, over a distance ~ Ljj, which is necessary 
for efficient acceleration of the bunch. With shorter focal 
length focusing and wq < A p , the interaction is char- 
acterized by filamented laser propagation and excessive 
curvature of the plasma wave, resulting in broad electron 
spectra. The experimental observations are supported by 
a series of 2D particle-in-cell (PIC) code simulations. 

The experiments were carried out on the Ti:Sapphirc 
Astra laser which provided pulses of energy up to 
E L = 700 mJ, pulse length r = 50(±5) fs FWHM at 
a central wavelength Ao = 795 nm. The laser was fo- 
cused with either f/3 or //12 off-axis parabolic mirrors, 
where the / refers to the ratio of focal length to initial di- 
ameter. The measured spot sizes (FWHM) at focus were 
4/im and 20 /im respectively, resulting in vacuum focused 
intensities of 5.5 x 10 19 Wcm~ 2 and 2.2 x 10 18 Wcm -2 
corresponding to ao of 5 and 1. The M 2 of both fo- 
cal spots was 2.6, giving zr (= ■kwq/M 2 \q) of 25 /xm 
and 600 /im respectively. The laser pulses were focused 
(z = 0) onto the front edge of a 3 mm diameter conical 
helium gas-jet. This provided initial electron densities of 



1 x 10 18 < 



< n e < 3 x 10 1 



The light transmitted through the plasma was mea- 
sured after reflection by a glass plate, which served to 
attenuate the intensity and limit spectral modifications 
caused by the window of the target chamber. Light 
within a //5 cone was collimated and re- imaged onto 
a 12 bit CCD camera. Light scattered at 90° to the di- 
rection of propagation was reimaged both perpendicular 
and parallel to the laser polarization plane. 

Electron beam divergence measurements were taken 
with a removable lanex scintillating screen with a 5° 
opening angle. The energy spectrum of the accelerated 
electrons was obtained with a magnetic spectrometer 
with a 5 x 1 mm steel collimator. The electrons were 
detected with image plates. These have a linear response 
over a dynamic range of 10 s , 50/mi spatial resolution, and 
allow calibrated, single-shot acquisition of the full spec- 
trum [15| . Alternatively, for the high charge monoener- 
getic beams, a lanex scintillating screen was used for 
high-repetition rate acquisition. 

With both f/3 and //12 focusing, at high density 
(cr > A p ) quasi-maxwellian electron energy spectra are 
observed. There is little shot-to-shot variation in spec- 
trum. This is indicative of self-modulation of the pulse 
fljj |. where the wakefield amplitude can reach the wave- 
breaking threshold over multiple wave periods. This, and 
the shorter dephasing length causes the large observed 
energy spread, as in Fig. QJi [4]. The divergence mea- 
surements show an almost uniform beam profile over an 
instrument limited opening angle of 5° as can be seen 
from Fig. [I]; i) and iii). Measurements taken with image 
plates closer to the interaction indicate that the actual 
divergence is significantly greater than this. 
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1: (color online) Electron spectra for a) f/3 at n e = 
3.0 x 10 19 cm" 3 and 2.4 x 10 19 cm" 3 ; b) f/12 at n e = 7.5 x 
10 18 cm -3 , c) Electron beam divergence measurements for i) 
f/3, n e = 3.0 x 10 19 cm -3 ii) f/3, n e = 1.5 x 10 1 " 
f/12, n e = 3.0 x 10 19 cm -3 iv) f/12, n e = 1.5 x 10 1 
The circles denote ~ 10 mrad intervals 



When cr < A p , the LWF regime is entered. In the case 
of the f/3 the number of accelerated electrons drops dra- 
matically, as can be seen in Fig.QJi ii). Coupled with the 
fall-off in electron flux due to the high beam divergence, 
the spectrometer can no longer detect any electron signal 
above noise level, even with the sensitive image plates. 

However, with //12 focusing, monoenergetic spectra 
are observed in the range 5.6 x 10 18 cm~ 3 < n c < 
2.0 x 10 19 cm~ 3 , Fig. QJ>. The electron beam divergence 
monitor consistently measures well collimated beams, 
with opening half- angle < 5 mrad, Fig. [Ik iv), in this mo- 
noenergetic regime. Large fluctuations in charge are mea- 
sured by the electron spectrometer. However, this can 
be explained by the electron beam pointing instability 
(3 x 10~ 3 sr) being greater than the collection solid-angle 
of the electron spectrometer collimator (5 x 10~ 6 sr). 

The variation in energy and energy spread may also 
be affected by this pointing instability. The ratio of the 
standard deviation to mean of the electron beam peak 
energy a/E = 0.18 for a set of 10 consecutive shots at an 
electron density of 7.5 x 10 18 cm -3 with E = 80 MeV. 

Side scattered images of the interaction indicate that 
the short focal length optic has a short channel length 
(~ zr) whilst for the long focal length they indicate ex- 
tended propagation. Fig.[2]shows images of the transmit- 
ted light at 800nm at z = 2 mm. For the f/3 (Fig. 2a) at 
all densities within the experimental range, the pulse is 
broken up into multiple filaments, even when the power 
for relativistic self-focusing — P cr — 17 ' .3{loq/lo p ) 2 GW- 
is exceeded, P/P cr > 1. As the density increased the 
average filament size decreased. 

The lack of guiding indicates why the f/3 (u> < A p ) 
cannot produce reproducible monoenergetic electrons un- 
der these conditions. To produce monoenergetic electron 
spectra the laser must propagate for distances on the or- 
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FIG. 2: (color online) Transmitted light at 800 nm imaged at 
z = 2 mm for (a) //3, n t = lx 10 19 cm" 3 and (d) f/12 at 
various densities. Also shown is the f/12 spot in vacuum (b) 
at focus and (c) z = 2 mm. (e) Mean spot width as a function 
of n e . Vertical line is P/P cr = 1, dotted line is wo = A p . 



der of the dephasing length. Due to the filamentation, it 
is unlikely that a single filament will be able to efficiently 
trap sufficient energy to generate beams of high charge 
and beam quality [171 ] compared with those obtained us- 
ing a long focal length optic. 

In contrast, using the f/12 optic (Fig. [2ji), the pulse 
emerges as a single filament in the density range of mo 
noenergetic electron spectra (6.7 — 13 x 10 18 cm -3 ). The 
spatial extent is comparable to that of the laser focus 
and scales with A p (Fig. 2e), indicating that the pulse 
has self-guided for a distance ~ 2 mm. Measurements 
of the transmitted energy suggest that at this point the 
pulse has depleted to below a sufficient power to maintain 
self-guiding. This is in contrast to the //3 case, where 
almost all of the energy is transmitted. 

For n e > 2 x 10 19 cm -3 , the spatial extent is large 
again. The depletion is so severe (> lOx) that self- 
guiding cannot be maintained over this length. The ob- 
served extended propagation allows the necessary non- 
linear pulse modification for injection of a single mo- 
noenergetic bunch of electrons. At the lowest densities 
(1.7 - 3.3 x 10 18 cm~ 3 ), for which P/P cr < 1, the spatial 
extent of the filament is close to that in vacuum at z = 2 
mm. This indicates that self-focusing only partially com- 
pensates for diffraction in this case. To our knowledge, 
this is the first evidence for self-guided propagation of 
short pulses (cr w \ p ). Note that previous results that 
suggested self-focusing is ineffective in this regime, only 
considered a short focal- length geometry (wq -C A p ) [l8[ . 

The two configurations were simulated using the 2D 
PIC code OSIRIS on an eight node Macintosh G5 cluster. 
The simulations were run with similar parameters to the 
experiments, for both the small and large focal spots, at 
4 particles/cell using a grid of up to 1500 x 6000 (> 18 
points/ A in x and up to 18 points/ A in y). A gaussian 



focus was used, due to the difficulty of initializing a pulse 
similar to the real focus. The embedded gaussian of the 
experimental focus was modeled to obtain similar diffrac- 
tion and Rayleigh lengths. This corresponds to 2 /im and 
10 fxm beam waists, defined as the half width at 1/e 2 
maximum at vacuum focus. Simulations with intermedi- 
ate focusing geometries were also performed. The density 
profile consisted of a 300 /im linear ramp followed by a 
plateau of density n e o to approximate the measured den- 
sity profile. The pulse started outside the plasma and its 
vacuum focus was set to be at the start of the plateau. 

Directly after focus, the 2 fim pulse develops a com- 
plex, strongly time-dependent mode structure. This elic- 
its a plasma response that enhances the structure via the 
self-focusing instability and evolves to a multi-filament 
structure (Fig. [3^) reminiscent of the analytical solution 
[l9| in the stationary approximation. Only about 30% of 
the pulse energy is trapped in a central filament. 

By contrast the 10 /^m pulse focuses radially and com- 
presses temporally, capturing almost all of the initial 
pulse energy (Fig. [3fc). The average pulse half- width 
shrinks to w 4 /zm, and the pulse propagates at this 
width for longer than than zr, until it is depleted. Un- 
like the 2 /Ltm case, the propagation is quasi-static and 
self-focusing occurs on a longer timescale. 

Electron density plots are also shown (Fig. [3ji, e) for 
both 2 /im and 10 /im simulations at the point where the 
plasma wave amplitude approaches its maximum value. 
Note that this amplitude maximum is at very different 
points in the simulation for the two pulses. Overlaid are 
contours of the laser eto, showing how the electron density 
profile directs the laser energy. 

The 2 ^m pulse interaction is shown close to focus af- 
ter which it breaks up and the wake amplitude decreases. 
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FIG. 3: (color online) 2D PIC simulations, showing: Contour 
plots of the laser pulse envelope (eE/mcujo) 2 for a) 2 /im, b) 
6 /im and c) 10 fim waist pulses after propagating 2zr. n e 
plotted for d) wo = 2 /j,m pulse at z = 0, e) »o = 10 fim 
pulse at z — 3zr. Overlaid are contours of the pulse at a = 
1, 1.5, 2, 3, 4. E L = 500 mj, r = 40 fs, n e = 2 x 10 19 cm -3 . 
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The electron trajectories from the edge of the pulse are 
attracted axially by the ion cavity produced by the in- 
tense pulse. The distance to reach axis should be similar 
to the distance the electrons slip back in the speed of 
light frame in this time, as they are rclativistic. However 
the pulse length in the case of the f/3 is much longer 
than the pulse width, and so the electrons move into the 
pulse envelope before being repelled by the ponderomo- 
tive force of the laser, as can be seen in Fig. [3jl. It is this 
crossing of the laser envelope by electrons that causes the 
break-up of the laser pulse. 

The 10 fiin pulse produces an initially more gently 
curved wake and reaches maximum amplitude later after 
non-linear modulation has increased its intensity. Even 
when the intensity is raised so that electron cavitation 
occurs quickly, there is still no breakup of the pulse as 
the width is similar to the pulse length. 

A scan in waist (between 1 and 14 fim) was simulated 
at n e0 /n c = 0.01, corresponding to 0.14 < w /X p < 2, for 
a fixed equivalent power, P = 13 TW. As noted above, 
larger (> 6 /jm) focal spots show stable quasi-static prop- 
agation (Fig. [3b, c), while spot sizes < 4 /im. (Fig. [3h) 
show beam break-up and dynamical self-focusing. Inter- 
mediate focusing geometries (~ 4 fim) show evidence of 
a transition region between the two regimes. 

An intensity scan (Jo/4 ^ I ^ 4Jo) at fixed density 
and a density scan 0.005 ^ n e /n c ^ 0.02 at fixed inten- 
sity were also simulated for both 2 fim and 6 /im geome- 
tries. The simulations were qualitatively invariant for 
both focusing geometries in this parameter range. The 
simulations also show that the pulse tends to self-focus 
to a diameter 2wo ~ X p and oscillate around an equilib- 
rium position (Fig. 3]). This is independent of the initial 
/ number, provided wq > X p . 

Despite the short pulse duration a simple stationary 
envelope treatment [201 ] can be used. For ao S> 1 the spot 
size w obeys d 2 w/d 2 t = V (l/w 3 - 16V2P/P cr ) for a lin- 
early polarized laser. Here Vq = [(cAo/7r)/(wQaQ)] 2 , W — 
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FIG. 4: (color online) Simulated temporal evolution of spot 
diameter at 1/e 2 (2w), with various focusing geometries for 
n e o/n c = 0.01, P/Pcr = 2. Vacuum foci are at z = /xm. 
Solid lines indicate simulated propagation. Dashed lines show 
theoretical propagation for a gaussian pulse in vacuum. 



w/woclq. Self- focusing will therefore start to dominate 
diffraction (d 2 w/dt 2 = 0) for w/X p = {\P/ P cr ) 1/e /vr. 
Hence for a mildly relativistic pulse, e.g. P/P cr « 2, 
the matched spot size 2w as (2/3) X p . This is in rea- 
sonable agreement with the simulations, the discrepancy 
being due to the assumptions of constant density and 
the do S> 1 approximation. For tight focusing geometries 
and high intensities where wo < X v and for ciq 3> 1, com- 
plete electron blowout can occur \l4 and a non-adiabatic 
kinetic analysis of self-focusing is required. 

In conclusion, we have demonstrated the self-guiding 
of an ultra-short laser pulse (ct ~ A p ) provided the 
vacuum spot-size wo > X p . Pulses with vacuum spot- 
size greater than this can propagate many zr, oscillating 
about a matched spot size. However, focusing to smaller 
than ~ A p will result in mode structures, which through 
sclf-focusing seeded beam break-up prevent production of 
high quality electron beams. So for reproducible produc- 
tion of monoenergetic electron, the ratio of pulse length 
to spot size should be < 1. To overcome low laser pow- 
ers, increasing pulse intensity through pulse compression 

is more 
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and photon deceleration of the leading edge 
effective than focusing more tightly. Careful choice of 
spot size, striking a balance between pump depletion with 
large / numbers and filamentation at low / numbers, will 
increase the efficacy of the interaction. This competition 
between over and under-focusing will be important even 
at higher intensities and so necessitates a careful choice 
of focusing conditions in future experiments. 
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